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Abstract 

A programming tactic involving polyhedra is reported that has been widely applied in the 
polyhedral analysis of (constraint) logic programs. The method enables the computations 
^— ^ ' of convex hulls that are required for polyhedral analysis to be coded with linear constraint 

, solving machinery that is available in many Prolog systems. 
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1 Introduction 

C/3 . 

O ' Polyhedra have been widely applied in program analysis IjCousot and Halbwachs 1978|l 

particularly for reasoning about logic and constraint logic programs. In this context 
polyhedra have been used in binding-time analysis ( |Vanhoof an d Bruynooghc 2001') , 
cdr-coded list analysis ( |Horspool 1990| ), argument-size analysis ( .Benoy and King 1996| ), 
5^ \ time-complexity analysis ( |King et al. 1997| ), high-precision groundness analysis IjCodish et al. 200Hl . 

type analysis ( |Saglam and Gallagher 1997| ), termination checking IjCodish and Taboch 1999|l 
and termination inference IjMesnard and Neumerkel 2001IIGenaim and Codish 2001|l . 

All these techniques use polyhedra to describe relevant properties of the pro- 
gram and manipulate polyhedra using operations that include projection, emptiness 
checking, inclusion testing for polyhedra, intersection of polyhedra (meet) and the 
convex hull (join). The classic approach to polyhedral analysis IjCousot and Halbwachs 1978|l 
uses two representations: (i) frames and rays and (ii) systems of (non-strict) lin- 
ear inequalities and employs the Chernikova algorithm to convert between them 
| |Le Verge 1992| |. The rationale for this dual representation is that the convex hull 
can be computed straightforwardly with frames and rays whereas intersection is 
more simply computed over systems of linear inequalities. A simpler tactic that has 
been widely adopted in the analysis of logic programs is to use only the linear in- 
equality representation and compute the convex hull by adapting | |Benoy and King 1996| ) 
a relaxation technique proposed in | |De Backer and Beringer 1993| |. The elegance 
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of this approach is that it enables the convex hull to be computed without re- 
course to a dual representation: the problem is recast as a projection problem that 
can be subcontracted to standard linear constraint solving machinery with mini- 
mal coding effort. Moreover, the performance is acceptable for many applications. 
In fact this technique has been widely applied in the analysis of logic programs 
HCodish and Taboch lOMIlGenaim and Codish"200Tl|King et al. 1997|IMesnard and Neumerkel 20011 
[Saglam and Gallagher 1997| ). The next section outlines the method and the follow- 
ing section, an example implementation. The final section presents the concluding 
discussion. 



2 Method 

Consider two arbitrary polyhedra, Pi and P2, represented in standard form: 

Pi = {xe Q" I Aix < Bi} P2 = {xe Q" | A2X < B2} 

such that Pi / and P2 7^ so that the problem is non-trivial. Note that AiX < Bi 
are non-strict and therefore Pi and P2 are both closed. The problem in essence is to 
compute the smallest polyhedron that includes Pi and P2 . Interestingly, the convex 
hull of Pi U P2 is not necessarily closed as is illustrated in the following example. 

Example 2.1 

Consider the 2-dimensional polyhedra Pi and P2 defined by: 



Pi = < f G 
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Observe that Pi ~ {(0, 1)} is a point whereas P2 = {{x,y) £ | x = y A < x } is 
a half-line. Note too that Pi and P2 are closed whereas the convex hull of Pi U P2 
excludes the points {(a;,y)GQ^|a;>OAy = a; + l} and hence is not closed (see 
the diagram below). 

Vi 

3- 






12 3 
Pi and P2 



0123a; 
convex hull of Pi U P2 



i 2 3 X 
closure of convex hull of Pi U P2 



Since the convex hull of Pi U P2 is not necessarily closed, the convex hull cannot 
always be represented by a system of non-strict linear inequalities; in order to 
overcome this problem, the closure of the convex hull of Pi U P2 is computed. The 
starting point for our construction is the convex hull of Pi U P2 that is given by: 



Ph = {x^ 



X — (JlXl + <T2X2 A 

Aixi < Bi A 



CTl + 0-2 = 1 
A2X2 < B2 



A 
A 



< (Ti A 

< (72 
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To avoid the non-linearity x — aiXi+a2X2 the system can be reformulated (relaxed) 
by putting yi — aiXi and y2 — (72X2 so that x = yi + y2 and Aii/i < CFiBi to define: 



PcH ^ {x£ 



X = yi+y2 A (Ti + (T2 = 1 A < CTi A 

^lyi < CTiSi A ^2^2 < CT2B2 A < (72 

Observe that Ph ^ Pch- Moreover, unlike Ph, Pch is expressed in terms of a 
system of linear inequalities. Note too that Pch is closed since the projection of a 
system of non-strict linear inequalities is closed. In fact the following proposition 
asserts that Pch coincides with the closure of the convex hull of Pi U P2 ■ 

Proposition 2.1 

Pch is the closure of the convex hull of Pi and P2 . 

The proof uses the concept of a recession cone. The recession cone of a polyhedron 
P, denoted 0+P, is defined by: 0+P = {y £ Q" | VA > . Vx e P . f -I- Ay G P}. 
The intuition is that 0+P includes a vector y whenever P includes all the half-lines 
in the direction of y that start in P. 

Proof 

Suppose P, = {f e Q" I Aix < B,}. Theorem 19.6 of l|Rockafellar 1970|l states that 
the closure of the convex hull of Pi U P2 is the set (0+Pi + P2) U (Pi + O+P2) U 
(U{(7iPi -I- (72 P2 I (71 -I- 0-2 = 1 A < (71, (72}). Intuitively, 0+Pi + P2 is P2 extended 
in the directions of half-lines contained within Pi. Let x ^ Pi, then y e O+P^ if 
and only if Ai{x + Xy) < Bi for all A > which holds if and only ii Aiy < 
URockafellar 1970 )[pp 62]. Therefore 0+Pi P2 e Q" | ^ = yi + ^2 A Aiyi < 

A ^2^2 < B2} and similarly Pi-hO+P2 = £ Q" | f = yi + 1/2 A Aiyi < 
Bi A A2y2 < 0}. Furthermore, U{(7iPi + (72P2 | (7i+(72 = 1 A < (71,(72} = 
{x G Q" I ai+(T2 = 1 A < (71,(72 A X yi +y2 A Aiyi < uiBi A ^22/2 < 0-2 P2}- 
Observe that {x e Q" \ x = yi+y2/\ Aiyi < aiBi A ^21/2 < (J2B2} coincides with 
the sets (i) O+P1+P2, (u) P1-I-O+P2 and (in) U{(7iPi-|-(72P2 | (7i-(-(72 = 1 A < (7i, (72} 
when (i) ui — Q and (72 = 1, (ii) ui ~ \ and (72 = and (iii) (7i+(72 — 1 and < (7i, (72 
respectively. Therefore Pch is the closure of the convex hull. □ 

This result leads to an algorithm for computing the closure of the convex hull: 
construct the systems Aiyi ^ <yiBi by scaling the constant vectors Bi by (7,, add 
the constraints x = yi -|- ^2, (7i -I- (72 = 1 and < cr^, then eliminate variables 
other than x using projection to obtain Pch in terms of x. Hence the closure of the 
convex hull can be computed without recourse to another representation. This is 
illustrated below. 

Example 2.2 

Returning to example 12. II consider the systems AiX < Bf. 



Pi ^\{x.y) e 



x<OA-x<0 A 
y < lA-y < -1 



I P2 = l{x,y)e' 



1,2 



X - y < A ' 
-X + y < A 
-X < 
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Adding a? = yi + 2/2, ci + (72 = 1 and < (Ti leads to the following system: 
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Eliminating the variables Xi, yi and leads to the solution: 

PcH = {{x,y) &Q^\0<xAx<yAy<x + l} 

Theorem 19.6 of IjRockafellar 1970(1 . which is used in the proof, asserts that Pch 
includes Pi + O+P2 = Pi + P2 = {{x,y) G Q"^ \ x > A y ^ x + 1} and therefore 
includes the points {{x,y) gQ^ \ x>OAy~x + l}, and hence ensures closure. 
Note that calculating Pch without the inequalities < (Ti and < (T2 - the 
relaxation advocated in ( |De Backer and Beringer 1993| ) for computing convex hull 
- gives {{x,y) G | < x} which is incorrect. 

3 Implementation 

This section shows how closure of the convex hull can be implemented elegantly 
using a linear solver in particular the CLP(Q) library IjHolzbaur 1995|l . The be- 
haviour of a predicate is described with the aid of modes, that is, + indicates an 
argument that should be instantiated to a non-variable term when the predicate is 
called; - indicates an argument that should be uninstantiated; and ? indicates an 
argument that may or may not be instantiated l|Deransart et al. 1996|l . 

3.1 Closed Polyhedra 

Closed polyhedra will be represented by lists (conjunctions) of linear constraints 
of the form c ::— e < e \ e — e \ e > e where expressions take the form e ::= x \ 
n \ n * X \ —e \ e + e \ e — e and n is a rational number and a; is a variable. A 
convenient representation for a closed polyhedron is a (non-ground) list of con- 
straints. This representation is interpreted with respect to a totally ordered (finite) 
set of variables. The ordering governs the mapping of each variable to its specific 
dimension. In practise, the ordering on variables is itself represented by the po- 
sition of each variable within a list. Specifically, if C is a list of linear constraints 
[ci, . . . , Cm] and X is a list of variables [xi, . . . , a:„], then the represented polyhedron 
is Pc,x = {(yi,...,y„) e Q" | (A^^^x, = y,) ho (A^LiCj)}. Note that although 
the order of variables in X is significant, the order of the constraints in C is not. 
Finally, let vars{o) denote the set of variables occurring in the syntactic object o. 

Example 3.1 

The polyhedron Pi from example 12.21 can be represented by the lists Ci = [x = 0, 
y = I] and X — [x,y\, that is, Pi = Pci,x- Moreover, P2 — Pc-2,x where 
C2 = [x — y, X > 0] or alternatively C2 = [y + z > x,x > y + 2 z,y > z > Q]. 
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Hence the dimension of Pc,x is defined by the length of the list X rather than the 
number of variables in C. 

3.2 Projection 

Projection is central to computing the convex hull. The desire, therefore, is to 
construct a predicate project(+Xs,+Cxs,-ProjectCxs) that is true when for a 
given list of dimensions Xs and a given list of constraints Cxs, ProjectCxs is the 
projection of Cxs onto Xs. The specification of such a predicate is given below. 

precondit ions : 

• Xs is a closed list with distinct variables as elements, 

• Cxs is a closed list of linear constraints, 

• Cxs is satisfiable. 

postconditions: 

• Xs is a closed list with distinct variables as elements, 

• ProjectCxs is a closed list of linear constraints, 

• wars(ProjectCxs) C i;ars(Xs), 

• -Pcxs.Xs = -PprojectCxs.Xs- 

Such a predicate can be constructed by adding the given constraints to the store 
and then invoking the projection facility provided in the CLP(Q) library, that is, 
the predicate dump (+Target , -NewVars, -CodedAnswer) (jHolzbaur 1995|l . Quot- 
ing from the manual: "[dump] reflects the constraints on the target variables into 
a term, where Target and NewVars are lists of variables of equal length and 
CodedAnswer is the term representation of the projection of constraints onto the 
target variables where the target variables are replaced by the corresponding vari- 
ables from NewVars". This leads to the following implementation of project: 

:- use_module (library (clpq) ) . 

project (Xs, Cxs, ProjectCxs) :- 
tell_cs (Cxs) , 

dump(Xs, Vs , ProjectCxs), Xs = Vs. 
tell_cs( [] ) . 

tell_cs([C|Cs]) :- {C}, tell_cs(Cs). 

Example 3.2 

For example, the query project ( [X, Z] , [X < Y, Y < Z] , ProjectCs) will cor- 
rectly bind Cs to [X-Z<0] . However, correctness of this predicate is compromised by 
existing constraints in the store. For instance, the compound query {X = Z + 1}, 
project ([X, Z] , [X < Y, Y < Z] , ProjectCs) will fail because constraints 
posted within tell_cs interact with those already in the store. 
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To insulate the constraints posted in tell_cs, both the variables Xs and the con- 
straints Cxs need to be renamed. Renaming is trivial with the builtin copy.term 
but care must be taken to ensure that Xs and Cxs are renamed consistently, that 
is that variable sharing in Xs and Cxs is preserved in the copies. However, in SIC- 
Stus Prolog copy_terin(Terin, Cpy) copies any constraints in the store that in- 
volve variables in Term. For example, the query {X=Y}, copy_terin(X=Y+l , Cpy) 
will bind Cpy to _A=_B+1 where _A and _B arc fresh variables. It will also copy 
the constraint X = Y by posting the new constraint Jl = _B to the store. To nul- 
lify this effect, copy_terin is called within the scope of calljresidue. The call 
call_residue(copy_term(X=Y+l , Cpy), Residue) residuates any new constraint 
into Residue instead of posting it to the store, thereby copying the term without 
copying any constraint. Whether residuation is required depends on the particular 
Prolog system. This leads to the following (SICStus Prolog specific) revision: 

project (Xs, Cxs, ProjectCxs) :- 

call_residue(copy_term(Xs-Cxs, CpyXs-CpyCxs) , _) , 

tell_cs (CpyCxs) , 

dumpCCpyXs, Vs, ProjectCxs), Xs = Vs. 
Example 3.3 

Using this revision, the query {X = Z + 1}, project ([X, Z] , [X < Y, Y < Z] , 
ProjectCs) will succeed binding Pro jectCs to [X-Z<0] . However, adding Z = 5 to 
the list of constraints induces an error. The problem is that posting the constraints 
binds Z to 5 so that dump is called with its first argument instantiated to a list that 
contains a non-variable term. 

A pre-processing predicate prepare_dump is therefore introduced to ensure that 
dump is called correctly. The following revision to project, in effect, extends the 
facility provided by dump to capture constraints over both uninstantiated and in- 
stantiated variables: 

project (Xs, Cxs, ProjectCxs) :- 

call_residue(copy_term(Xs-Cxs, CpyXs-CpyCxs) , _) , 
tell_cs (CpyCxs) , 

prepare_dump(CpyXs, Xs, Zs, DumpCxs, ProjectCxs), 
dump(Zs, Vs, DumpCxs), Xs = Vs. 

prepare_dump( [] , [] , [], Cs, Cs) . 
prepare_dump ( [X I Xs] , Ysin, ZsOut, Csin, CsOut) :- 
(ground (X) -> 

YsIn = [YlYs] , 

ZsOut = [_ I Zs] , 

CsOut = [Y=X|Cs] 

YsIn = LlYs] , 
ZsOut = [X I Zs] , 
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CsOut = Cs 

), 

prepare_dump(Xs, Ys , Zs, Csin, Cs) . 

The literal prepare_dump(+Xs, +Ys, -Zs, ?CsIn, -CsDut) is true for a given list 
Xs which contains either variables or numbers (or a mixture of the two) and a given 
list Ys which contains only variables, if 

• Zs is the list obtained by substituting the non- variable terms of Xs with fresh 
variables and 

• CsOut is an open ended list of equality constraints with CsIn at its end 
that contains one equality constraint for each number in Xs. Each constraint 
equates a numeric element of Xs with the element of Ys that is in the same 
list position. 

Thecanprepare_dump([Xl, 1, X3, 2], [A, B, C, D] , Zs, CsIn, CsOut), for 
instance, wiU bind Zs to [X1,_A,X3,_B] and CsOut to [B=l ,D=2 I CsIn] . The pred- 
icate ensures that dump is called with its first argument bound to a list of free 
variables even when the list Xs includes numbers. In the CLP(Q) library, numbers 
coincide with rationals which are represented as compound (ground) terms of the 
form rat(n, d) where n and d are integers. The ground (X) test effectively checks 
whether X is instantiated to a number; the test number (X) is inappropriate since it 
would always fail. 

Example 3.4 

Consider again example 13. II The second representation of P2 can be simplified by 
using projection as follows: 

I ?- Cs = [Y+Z>=X,X>=Y+2*Z,Y>=0,Z>=0] , proj ect ( [X , Y] , Cs, ProjectCs) . 

ProjectCs = [Y>=0,X=Y] ? ; 

no 

The system Cs is expressed over 3 variables and therefore defines a 3 dimensional 
space. Intuitively, the projection onto [X, Y] is the shadow cast by Pcs,[x,y,z] onto 
the 2 dimensional space over X and Y. The projection ProjectCs in fact defines 
a half-line confined to the first quadrant since, by rearranging Cs, it follows that 
Pcs,[x,Y,z] = {{x, y, z) e I X = y A < y A z = 0}. 

3.3 Convex Hull 

The specification for the main predicate convexJiullC+Xs , +Cxs, +Ys, +Cys, 
-Zs, -Czs), and then its code, is given below. 

preconditions: 

• Xs is a closed list with distinct variables as elements and likewise for Ys, 

• Xs and Ys have the same length, 

• varsiXs) fl vars{Ys) — 0, 

• Cxs and Cys are closed lists of linear constraints, 

• Cxs and Cys are both satisfiable, 

• varsiCxs) C vars(Xs) and vars{Cys) C vars(Ys). 
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postconditions: 

• Xs, Ys and Zs are closed lists with distinct variables as elements, 

• Zs is the same length as both Xs and Ys, 

• Czs is a closed list of linear constraints, 

• vars{Czs) C vars(Zs) and (uars(Xs) U vars{Ys)) r]vars{Zs) = 0, 

• -Pczs.zs is the closure of the convex hull of Pcxs.xs U Pcys,Ys- 

coiivex_hull (Xs , Cxs, Ys, Cys, Zs, Czs) :- 
scale (Cxs, Sigl, [] , Cls) , 
scaleCCys, Sig2, Cls, C2s) , 
add_vect (Xs , Ys, Zs, C2s, C3s) , 

project (Zs, [Sigl >= 0, Sig2 >= 0, Sigl+Sig2 = l|C3s], Czs). 

scale ( [] , _, Cs, Cs) . 
scale([Cl|Cls] , Sig, C2s, C3s) :- 

CI =. . [RelDp, Al, Bl] , 

C2 =. . [RelOp, A2, B2] , 

mul_exp(Al, Sig, A2) , 

inul_exp(Bl, Sig, B2) , 

scale(Cls, Sig, [C2|C2s], C3s) . 

mul_exp(El, Sigma, E2) :- once(mulexp(El , Sigma, E2)). 



mulexp ( X , _ , 


X) 


- var(X) . 






mulexp(N*X, _, 


N*X) 


- ground (N) , var(X) . 






mulexp ( -X, Sig, 


-Y) 


- mulexp (X, Sig, Y) . 






mulexp (A+B, Sig, 


C+D) 


- mulexp(A, Sig, C) , mulexp(B, 


Sig, 


D) 


mulexp (A-B, Sig, 


C-D) 


- mulexp(A, Sig, C) , mulexp(B, 


Sig, 


D) 


mulexp ( N, Sig, 


N*Sig) 


- ground (N) . 






add_vect ( [] , [] , 


[], Cs, 


Cs) . 






add_vect( [UlUs] , 


[ViVs] , 


[WlWs] , Cls, C2s) :- 






add_vect (Us , 


Vs, Ws, 


[W = U+VlCls] , C2s) . 







The predicate mulexp (?E1 , ?Sigma , -E2) scales the numeric constants that oc- 
cur within El by the variable Sigma, providing they are not coefficients of variables, 
to obtain the expression E2. Note that Sigma is a variable and the expression El may 
be a variable, hence both El and Sigma have mode ? rather than +. Since a non- 
ground representation is employed for expressions, the test var (X) is used to deter- 
mine whether the expression is a variable. As before, the test ground (N) detects nu- 
meric constants - rational numbers - which are the only type of subexpressions that 
are ground. Observe that mulexp can return more than one solution, for example, 
mulexp(X, Sig, E2) generates E2 = X; X = -(_A), E2 = -(_A); X = -(-(_A)), 
E2 = -(-(JV)) etc as solutions. Thus the pruning operator once is applied within 
mul_exp(?El, ?Sigma, -E2) to prevent erroneous solutions. 
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The predicate scale (+Cls , ?Sigma, ?C2s, -C3s) scales each constraint within 
the hst Cls by the variable Sigma. Each constraint consists of a binary operator and 
two expressions, and scaling is applied to the numeric constants in each expression as 
specified by mul_exp. For example, scale ( [X+2 >= 1+Y, Y = Z] , Sigma, Tail, 
ScaledCs) binds ScaledCs to [Y = Z, X+2*Sigma >= l*Sigma+Y I Tail] . Note 
that scale finesses the problem of putting Cxs and Cys into the standard form 
■^iVi < Bi before applying scaling. In standard form, X+2 >= 1+Y is Y-X =< 1 but 
scaling constants on both sides of the relational operator preserves equivalence in 
that X+2*Sig >= l*Sig+Y is equivalent to Y-X =< l*Sig. The use of a difference 
list avoids an unnecessary call to append in the body of convexJiull. 

The predicate add_vect (+Us , +Vs, -Ws, ?Cls, -C2s) operates on the lists Us 
= [Ui , . . . , U„] and Vs = [Vi , . . . , V„] which correspond to the vectors yi and 
?/2 (as introduced in section ^ . The argument Ws is instantiated to another list 
of variables [Wi , . . . , W„] , which corresponds with x. The predicate creates the 
system of equalities [Wi = Ui+Vi, W„ = U„+V„] corresponding to the system 

X = yi + 2/2- The scaled constraints output by the two calls to scale are passed to 
add_vect via its accumulator and thereby combined with the system of equalities. 
For example, the caU add_vect( [XI, Yl] , [X2, Y2] , Ws, Tail, Cs) returns the 
bindings Cs = [_A=Y1+Y2,_B=X1+X2 iTail] and Ws = [_B,_A]. 

The predicate convexJiull (Xs , Cxs, Ys, Cys, Zs, Czs) takes, as input, two 
lists of constraints (Cxs and Cys) and their corresponding lists of variables (Xs 
and Ys) and produces as output a single list of constraints Czs over the vari- 
ables Zs that represents the closure of the convex hull of the two input poly- 
hedra. If Xs and Ys are not variable disjoint, then the pre-requisite can be sat- 
isfied by appropriately renaming variables. Specifically, the variables Xs and con- 
straints Cxs can be renamed with copy_term(Xs-Cxs , CpyXs-CpyCxs) and the call 
convexJiull (Xs , Cxs, Ys, Cys, Zs, Czs) replaced with convex Jiull (CpyXs , 
CpyCxs, Ys , Cys, Zs, Czs). Since the integrity of the constraint store is pre- 
served by project and since project is the only source of interaction with the 
store, then it follows that convexJiull also does not side-effect any existing con- 
straints. The following is an illustrative example. 

Example 3.5 

Running this code on the data of Example 12 . 21 gives: 

I ?- convex_hull([Xl,Yl] , [X1=0,Y1=1] , [X2,Y2] , [X2>=0, Y2=X2] ,V,S) . 
S = [_A>=0,_A-_B>=-1,_A-_B=<0] , 
V = [_A,_B] ? ; 
no 

4 Discussion 

This section discusses the method proposed in the paper, comparing it with related 
techniques. The Chernikova method is exponential in the worst-case ( |Le Verge 1992| | 
and the Fourier-Motzkin method, like all projection techniques over linear inequal- 
ities l|Chandru et al. 20n0|l . is also exponential. The exponential behaviour of both 




Fig. 1. (i) Pi and P2, (ii) Qi and Q2, (iii) conv{Pi U P2), (iv) conv{Qi U Q2) 

methods stems from the same source: the possibly exponential relationship between 
the number of vertices and the number of half-spaces that define a polyhedron. In 
fact the problem of calculating the closure of the convex hull of two polyhedra is also 
exponential even for bounded polyhedra (polytopes). This can be demonstrated by 
considering the so-called cross polytope in n-dimensions which is the polyhedron 
with the vertex set {(±1, 0, . . . , 0), (0, ±1, . . . , 0), . . . , (0, 0, . . . , ±1)}. The cross poly- 
tope can be defined by no less than 2" inequalities yet can arise as the convex hull 
of two polyhedra both of which can be defined with 0{n) inequalities. Specifically 
consider the n-dimensional polyhedra 

Pi = {(Xi, . . . , Xn) e Q" I (Etl < 1) A {/\]=lX, < 0)} 

P2 - {(xi, . . . , Xn) e Q" I (Er=i < 1) A (Aj'^i - xj < 0)} 

Because Pi and P2 are polytopes, they can be expressed in terms of their vertices: 

Pi =con«({(0,0,...,0), (-1,0,..., 0), (0,-l,...,0), (0,0,...,-l)}) 
P2 = con«({(0,0,...,0), (1,0,..., 0), (0,1,...,0), (0,0,...,1)}) 

Since (0, 0, ... , 0) is convexly spanned by (1,0,..., 0) and (—1, 0, . . . , 0), it follows 
that d{conv{Pi U P2)) = conv(Pi U P2) = ^^^({(±1, 0, . . . , 0), (0, ±1, . . . , 0), . . ., 
(0, 0, ... , ±1)}) which is the n-dimensional cross polytope. The 2 and 3 dimensional 
cases are denoted in Figure ^ by (i) Pi and P2 and (ii) Qi and Q2 respectively for 
which the cross polytopes are a solid square and an octahedron. Hence the problem 
of calculating the closure of the convex hull is intrinsically exponential irrespective 
of the algorithm employed. 



Example 4-1 



Computing Convex Hulls with a Linear Solver 
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The following query illustrates how the hull algorithm yields an exponential number 
of inequalities for the 4 dimensional case. 



I ?- Xs = [XI, X2, X3, X4] , Ys = [Yl, Y2, Y3, Y4] , 

Cxs = [-1 =< X1+X2+X3+X4, XI =< 0, X2 =< 0, X3 =< 0, X4 =< 0] , 
Cys = [ Y1+Y2+Y3+Y4 =< 1, =< Yl, =< Y2, =< Y3, =< Y4] , 
convex_hull (Xs , Cxs, Ys , Cys, Zs, Czs) , 
Zs = [A, B, C, D] . 

Czs = [A-B+C+D>=-1, A+B-C-D=<1, A+B+C+D>=-1, A-B-C-D=<1, 

A-B-C+D>=-1, A+B+C-D=<1, A+B-C+D>=-1, A-B+C-D=<1, 

A-B+C-D>=-1, A+B-C+D=<1, A+B+C-D>=-1, A-B-C+D=<1, 

A-B-C-D>=-1, A+B+C+D=<1, A+B-C-D>=-1, A-B+C+D=<1] ? ; 



However, it would be wrong to conclude from these examples that the frame and ray 
representation is preferable - inequalities are unavoidable since they are required 
for other polyhedral operations. 

Despite the scaling problems that are inherent to any convex hull algorithm, in 
practise the technique proposed in this paper has been widely applied in logic pro- 
gramming (jCodish and Taboch 1999: Genaim and Codi sh 20011 [King et al. T997| 
IMesnard and Neumerkel 2 001 Saglam and Gallagher 1997| ), mostly to satisfaction. 
For example, in the context of inferring termination conditions for logic programs 
this method is feasible since it accounts for 42% of this first pass of the analysis and 
the first pass itself constitutes only 23% of the total analysis time (jMesnard and Neumerkel 200 l|l . 
Whether the approach presented in this paper is applicable depends on the appli- 
cation context. When only standard domain operations are required and perfor- 
mance is not an issue, this method has much to commend it. However, when the 
application has to additionally reason, say, about integral points IjAncourt 19911 
IQuinton et al. 1997| ) or parameterised polyhedra IjLoechner and Wilde 1997|l then 
specialised polyhedral libraries are required. Further, if performance is important, 
then recourse should be made to a polyhedral library, since a state-of-the-art im- 
plementation employing the Chernikova algorithm ( |Bagnara et al. 2002| ), will out- 
perform the approach presented here. 

We have presented a Prolog program for computing convex hulls using linear 
solver machinery. As Holzbaur's library is also available for CIAO Prolog, ECLiPSe, 
XSB and Yap Prolog, the technique can be easily adapted to these systems. The 
method is a reasonable compromise between conciseness, clarity and efficiency and 
variants of this program have now been widely deployed. 
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